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Abstract. We compare two QCD-inspired quark models with four-fermion inter- 
action, without and with the remnant coupling to low-energy gluons, in the regime 
of dynamical chiral symmetry breaking (DCSB). The first one, the Nambu- Jona- 
Lasinio (NJL) model ensures the factorization of scalar and pseudoscalar meson poles 
in correlators, the well-known Nambu relation between the scalar meson mass and 
the dynamical quark mass, m a = 2m^ yn , and the residual chiral symmetry in cou- 
pling constants characteristic for the linear cr-model. The second one, the Gauged 
NJL model, happens to be qualitatively different from the NJL model, namely, the 
Nambu relation is not valid and the factorization of light meson poles does not entail 
the residual chiral symmetry, i.e. it does not result in a linear cr-model. The more 
complicated DCSB pattern in the GNJL model is fully explained in terms of excited 
meson states with the same quantum numbers. The asymptotic restrictions on pa- 
rameters of scalar and pseudoscalar meson states are derived from the requirement 
of chiral symmetry restoration at high energies. 
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1. Introduction 



The QCD-inspired quark models with attractive four-fermion interaction are often 
applied for the truncation of the low-energy QCD in the hadronization regime [1-9]. 
In the common approach [1,2] the local four-fermion interaction only is involved to 
induce the DCSB due to strong attraction in the scalar channel. As a consequence, the 
dynamical quark mass m dyn is created, massless pions (in the chiral limit, m current = 
0) and the massive scalar meson with the mass, m c = 2md yn arise and the residual 
chiral symmetry characteristic for the linear a-model [2] holds in coupling constants 
(see below). On the other hand the effective QCD action for quarks at low energies 
contains, in general, the remnant interaction to low-energy, background gluons which 
is responsible for the confinement and hadron properties at intermediate energies, in 
particular, for the formation of excited meson states [10]. Therefore, when gluons are 
treated as being low-energy, background ones the Gauged NJL model may be thought 
of as a better truncation of the low-energy QCD effective action [5]. As the DCSB 
holds also for such a model the appearance of large dynamical quark mass gives rise 
the possibility to calculate nonperturbative gluon corrections for hadron parameters 
in the condensate approach [3-6]. 

When comparing two models in the description of hadron properties one has to 
examine to what extent the Gauged NJL model leads to the same hadron lagrangian 
as the conventional NJL one, after proper renormalization of its basic characteristics 
[4]. It is one of the aims of our paper to show that two models are qualitatively 
different since the coupling to background gluons breaks the linear a-model relations 
for masses and coupling constants which are valid in the non-gauged case. 

In order to explain the roots of this difference we consider the qualitative features 
of QCD hadron spectrum in the planar limit (large N c ). In this approach the corre- 
lators of color-singlet quark currents are saturated by narrow meson resonances. In 
particular, the two-point correlators of scalar and pseudoscalar quark densities are 
represented by the sum of related meson poles, 

n S (Q) = -J d 4 X exp(iqx) < T (tjj^(x) >planar 

II P (Q) = J d 4 x exp(iqx) < T (ip^(x) ^75^ (0)) > p i an ar 

in the Euclidean momentum region, q 2 = —Q 2 < 0. The absence of physical thresh- 
olds for q 2 > in the large N c limit is ensured by the confinement of color intermediate 
states. The high-energy asymptotics is controlled [11] by the perturbation theory and 
the operator product expansion due to asymptotic freedom of QCD. Thereby the same 
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correlators have the power-like behavior at large Q 2 , 

N O 2 

n s (g)| Q2 _ ~ n p {Q)\ Q ^oo * ^Q 2 in^V (2) 



When comparing (1) and (2) one concludes that the infinite series of resonances 
with the same quantum numbers should exist in order to reproduce the perturbative 
asymptotics. Besides, one can prove [12] that 



(n»(0)-n'(0)) =o U , (3) 



and the chiral symmetry is restored at high energies. 

Thus the QCD-inspired quark models are expected to reproduce the part of QCD 
meson spectrum in the planar limit and the asymptotic chiral symmetry restoration 
for higher energies. In particular, the conventional NJL model describes one scalar 
and one pseudoscalar (multiplet of) state whereas the Gauged NJL model should 
contain additional excited meson states with the same quantum numbers. Just the 
appearance of additional poles in correlators makes two models to be different in the 
DCSB pattern. 

The rest of the paper contains the derivation of two-point correlators in the NJL 
and Gauged NJL models, the proof of the violation of linear a-model relations and 
the soft-momentum determination of scalar meson parameters. At the end the two- 
resonance ansatz is performed to saturate the deviation of NJL and GNJL models. 
The asymptotic sum rules based on the chiral symmetry restoration are obtained for 
the resonance characteristics. 

2. The definitions 

We start from the NJL model without gluons and introduce the auxiliary scalar, 
a = a(x) +m c iyn', ijndyn = m ) and pseudoscalar, ir(x) fields. Respectively the basic 
lagrangian density can be presented in two forms (the euclidean-space formulation is 
taken here), 

C(x)) = $(id + iS(x) + lb P(x)) i> + ^- [(^) 2 - $75^) 

$(D + iS(x) + ^P{x))i) + ^ (V(x) + m) 2 + 7t 2 (x)) , (4) 

where ip = ipi stands for color fermion fields with iV c components, S(x),P(x) are 
external scalar and pseudoscalar sources and the Dirac operator is 

D = id + im + icr(x) + 7r(x)7 5 ; 8 = 1^.8^ (5) 

The vector and axial- vector fields are not included since the main story is interplayed 
between scalar and pseudoscalar channels. We simplify our analysis setting up the 
number of flavours N F = 1 and the current quark mass m q = 0. 
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Let us make the change of fields, a ^ a — S; tt ^ n — P. Then the regularized 
fermion vacuum functional is described by the following relations, 

In Z^ eg (a, it) = In ^exp^— J d 4 x C F (a(x), n (x) )^ _ = Af c Tr In _D| re9 (6) 

The corresponding effective action for auxiliary fields reads 



S e ff = -A^ c TrlnL)| reg 
N, 

9 



N r 

+ J (fx {a 2 + 2ma + n 2 - 2aS - 2mS - 2nP + S 2 + P 2 + m 2 ) (7) 



where a momentum-cutoff regularization is implied. The parameter g is a four-fermion 
constant. When relating to the BRZ denotations [8] it happens to be g 2 = 8tt 2 G s /A 2 . 

3. Parameters of scalar and pseudoscalar mesons in the conventional NJL 
model 

As usual the dynamic mass is subject to the mass-gap equation in which A is a 
0(4)-invariant cutoff, 

/ . 9 87T 2 9 , A 2 , ^.m 2 
m(A — m In — -) + 0(— ) = 0. 8 

g 2 m 2 ' y JS? J w 

Herein and further on we omit the terms of order 1 /A 2 which do not change the result 
drastically. This equation provides the cancellation of the linear in a terms in (]?]), 
with the tadpole in the fermion determinant, 

2N c m — L iN, 



— = K#>= (vi) Tr K es '-» = 

= W c [ ^ ™ (9) 
J P <A 16n 2 p 2 + m 2 

The quadratic in fields part of S e fj determines the inverse propagators ("kinetic 
terms") for meson fields, 

Seff = l/0i [°(Q)r*(QX-Q) + AQ)^(Q>(-Q)} , (10) 

where in the simple NJL model without gluons one has, 



g 2 L y p <A(2vr) 4 y 2 J v 2 " 

N ° / 77rT7tn7 5 A(— )«7 5 A(— ), (11) 



g 2 t y p <A(2vr) 4 10 y 2> 10 v 2 
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and the fermion propagator is given by A(Q/2) = (p + Q/2 + iro) 1 . The correspon- 
dence to the BRZ denotations [8] is following, 

1 „ 1 



n, 



i 

9s 



N C A 2 



9s 
2Nc 
9 2 



r CT ; n P 



< qq > 



m 



9s_ 

.<1p1p >£ 



m 



(12) 



Both integrals in (|TT|) are divergent. If we take into account the mass gap equation 
then they are logarithmically divergent only. Let us evaluate the integrand in (p|), 



Q 



Q 



tr(l, ij 5 ) [p + — + im] (l,ry 5 ) p- — + 



im 



4 ( (p, + %) 2 + m s 



Qt 



2 Q - 2 
p — — =p m 



fa~^) + m 



-i 



/ QfJ.\2 i 2 



-1 



- 2 (Q 2 + 2m 2 ± 2m 2 ) (fa + ^) 2 + m 



-1 



(13) 



The sum of demoninators is, 2p 2 + 2m 2 + (Q 2 /2) and after its separation in the 
numerator of Eq. (|T3"|) one arrives to the required decomposition of loop integrals. 
Two first terms at zero-momentum represent the tadpole integrals in ([|) and therefore 
precisely cancel the first terms in (|TT|). This is a work of the mass-gap eq.([8|). Thus, 
depending on the regularization we derive the inverse cr-model propagators, 



IV = (4 m 2 + Q 2 )I(Q 2 )+ZQ 2 + 0(-) = (ml + Q 2 )I a (Q 2 ) 
T n = Q 2 l(Q 2 ) + ZQ 2 + (-^) = Q 2 I 7T (Q 2 ), 



(14) 



where the parameter £ is a finite constant characteristic to the regularization, for in- 
stance, in the symmetric cutoff regularization (|TTD one has £ = N c /32tt 2 . Respectively 
the decay formfactor is provided by the last term in (O), 



HQ 2 



2N r 



d A p 



1 



(15) 



l P <A (2vr) 4 (p + \Q) 2 + m 2 (p- \Q) 2 + m 2 ' 

It is regularization independent up to a redefinition of A when neglecting the orders 
of 1/A 2 . For a convenience we display the analytic representation for I(Q 2 ), 



HQ 2 



8tt 2 

No 
8tt 2 



' A 2 
In — 



1 



Q 2 



dx In 1 + x(l — x) — - 



In 



A 2 



+1+W1+ 



Am 2 



rrr 



Q 2 



In 



1 + 



1 + ^ + 1 



(16) 
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for euclidean momenta Q 2 > 0. If one defines the regularization where £ = (or an 
appropriate subtraction of £-term in the arbitrary regularization) then one arrives to 
the residual chiral symmetry and thereby to the generalized linear cr-model, 

h{Q 2 ) = i,{Q 2 ) = i{Q 2 )- (17) 

Then as a consequence of the NJL-mechanism of DSB we have precisely vn a = 2m. 
The generating functional for scalar and pseudoscalar correlators can be derived after 
one integrates over a, ir variables keeping in the effective action the quadratic terms 
only (the main large- N c contribution). The gaussian integral is calculated at the 
extremum, 

and results in the generating functional, W = — In Z(S, P), 



(2vr) 4 



k-g)p^r--i^U(g) 

2 \ m l m I 



■\p{-Q) f^^r; 1 - £±t>) p( Q ) - 5(0)(2 7 r 4 )^)(g ) z < ^ > 

2 \ m A m / 



.(19) 



The calculation of two-point correlators correspond to the second derivative of this 
action with subtraction of disconnected product of condensates (it corresponds to 
elimination of the last term ~ S(0) in above Eq.). 

d A x exp(iQx) < T(^jji(}(x)^^(0)) > — < ipip >< ipip > 

^ 4 xexp(zgx) ^ ( ^ (o) (W^ + J d A yS{y)i < ^ >) ■ 

It is saturated by first terms in ([19]) quadratic in S(x). One has the scalar-meson 
pole and another constant term which plays some role [8] when one reduces linear 
sigma model to the chiral one. Similarly we obtain the pseudoscalar correlator which 
contains also the pole and the constant term. One can compare them to the correlators 
(124)-(127) and (152) of the BRZ paper and find the one-to-one correspondence if to 
go back to the Minkowski space where —% < ipip >e^< >a/=< QQ >■ The chiral 
symmetry restoration at high energies is ensured by (ITJ 



4. Soft-momentum predictions 

Let us now develop the soft-momentum expansion of above kinetic terms, 

r° = al + Q 2 a\ + • • • , T° = Q 2 a\ + Q A a n 2 + ■■■. (20) 
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The coefficients aj'^ are related to the Taylor expansion of (|T 
a° = m%(0), al = I o (Q)+m 2 M0) 

aj = 4(0), al = I'M (21) 

We see that in the vicinity of Q 2 ~ the would-be scalar-meson decay constant a\ 
and pion constant a\ do not coincide. Moreover their difference represents the finite 
constant nearly independent on a regularization, a\ — a\ = 4m 2 /' (0). We remark 
that the corresponding integral I'(0) is convergent and in the limit A — > oo takes the 
value I'(0) = —N c /487r 2 m 2 . On the other hand, when the residual chiral symmetry 
( [TTD is present (for £ = 0) the latter one reveals itself in a series of cr-model identities 
between the coefficients of soft-momentum expansion, 

4 = ™*, aW = K) 2 + «>••■ (22) 

We emphasize that the first relation is meaningful iff the next ones are valid. But only 
if I a = I w one can use the first relation in eq. (|22|) for the definition of the scalar meson 
mass. Otherwise it is misleading. Obviously, in the conventional NJL model the scalar 
meson mass is easily determined from the soft-momentum (derivative) expansion. In 
the next section we shall see that it is not the case for the Gauged NJL model. 

5. Gauged NJL model and status of a-model 

The gauged NJL model is prepared conventionally with help of the long derivative 
in the Dirac operator, id^ — > = id^ + G^ and G M = g s T a G^ is a soft gluon field. 
The latter leads to the following modification of quark propagators in loop integrals, 
eq-©, 

A(±|)= UQ + id + G + im^ . (23) 

Respectively we examine the main large- iV c contribution into the effective action given 
by averaging of functionals r CTj7r (A) over gluons, 

IV =< rV(A) > G . (24) 

We assume that the main DSB mechanism is due to the strong four-fermion interac- 
tion and soft gluons bring the corrections accumulated into condensates. Thereby we 
imply that the derivative or large-mass expansion is reasonably good in calculating 
the gluon contributions. 

Let us analyze the general structure of Y a ^ which defines the inverse propagators 
for a and 7r-fields according to (]To|) . For this purpose we first evaluate the integrand 
in ©,(0), 

Tr(l,i T5 ) (Q+id + G + im\ (l,z 75 ) +id + G + imj = 
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x ({id» + G») 2 + F-^Tm 2 + \%D]} 

x + i9 + Gfi f + F + m 2 ) 1 (25) 



where F = (1/4) [7^, 7^] ■ [D^,D U ]. The sum of demoninators is, 

£ = 2(Z^) 2 + 2F + + 2m 2 (26) 

and after its separation in the second multiplier of product (^) one arrives to the 
following decomposition of loop integrals (compare with eq. (p"4"D ) : 

T a = (Q 2 + 4m 2 )/ G (Q 2 ) + Q 2 ~(Q 2 ) + 0(^), 

K = Q 2 (l G (Q 2 )+E(Q 2 ))+0(^), (27) 
where now the employed functions are defined as follows, 



-(tr <x\ l-^-i /o,, 



W) = ~(tr <x| (^ + ^ + G M ) 2 + F + m 2 



x f(p M -% + ^ M + G M ) 2 + F + m 2 ) ^ 



(Q 2 ) = e-^tr<x|((^ + ^ + G M ) 2 + F + m 2 ) 



q 2 ^ J V 2 
1 /■ rf 4 p . (. , g„ 



-1 



x^[g,P + ^](fe t -^ + ^ + G M ) 2 + F + m 2 ) ) G . (28) 

The integral for S is not divergent. We have used the modified gap equation and the 
1/A 2 expansion which results in the independence of parameter £ on gluon conden- 
sates ( £ = N c /32tt 2 in the sharp cutoff regularization) . The representation for Iq 
and S in Eqs.(p8[) contains the trace over spin and color indices. In the second, in- 
tegral representation of above functions the derivatives act on gluon fields when one 
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expands, for example, those functionals in the perturbation series. In other words, we 
have to calculate the mixed matrix element < x\ ■ ■ ■ \ k = > . Evidently the relations 
of the linear cr-model (|T7|) would survive if it were true that, 



(Q 2 + 4m 2 )/ G (Q 2 ) + Q 2 ~(Q 2 ) I (Q 2 + m 2 a ) (l G (Q 2 ) + ~(Q 

4m 2 — m 2 
m 



or ~(Q 2 ) 1 ""„ 2 "' g / G (g a )- (29) 



By comparing of integrals in ([28D one can convince oneself that for arbitrary momenta 
these functions do not coincide, the equation (|29| ) is not valid and therefore the linear 
cr-model relations ([17|), (|22"D are not fulfilled. Thus we cannot effectively use the soft- 
momentum expansion for the calculation of scalar meson mass and its decay constant. 
Now let us consider the soft-momentum and large-mass expansion which leads to the 
expansion in terms of gluon condensates. We stress that the coefficients of large- 
mass expansion containing gluon condensates are free of ultraviolet divergences and 
thereby they do not depend on regularization if one neglects terms 0(1 /A 2 ). The 
new mass-gap equation takes the form, 

t2 8tt 2 2l A 2 <G 2 > 

A 2 - — - m 2 In — + — f- = 0. 30 



First terms of above expansion for T a ,T n read (see eq . fl21|) ) in the momentum- cut off 
regularization after applying the mass-gap equation: 



«o 




4 = ; -i • (si) 

2 48vr 2 m 2 1 on v 1 

where < G 2 U >= g 2 QCD < (G* u ) 2 >, £ = 1/4 = £ • 8tt 2 /N c . One can verify 
explicitly that the linear cr-model relations (PTj) , (P0|) are not fulfilled even after the 
renormalization to the chirally symmetric point, % ,7r — > a° ,7T — £; £ = N c /32tt 2 . 
For the completeness, we display also the soft-momentum expansion for functions 
I G (Q),E(Q), 

t (o) = A2 i i t I <G l»> Q 2 ( l | < G %> \\. 

GW> 8n 2 { n m 2 5 24A^ c m 4 m 2 ^6 40N c m 4 / / ' 

1 <G 2 > Q 2 <G 2 > 
8vr 2 24m 4 8vr 2 m 2 60m 4 ' 1 ; 
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Thus one can see that the relations ( 29f) cannot be fulfilled for any reasonable choice 
of £ (which gives 1(0) ^ 0). 

Let us refer our analysis to the identities (199)-(201) in the paper [8] for one- 
fermion loop (OFL) polarization operators, scalar, Tl s , pseudoscalar, Tip, and longi- 
tudinal axial- vector one, r4 0) . They are valid in the non-gauged NJL model. The 
PCAC relations (199), (200) connect the OFL pseudoscalar correlator, Tip, with the 
OFL longitudinal axial- vector one, 5^. In the GNJL model only these two rela- 
tions survive. Meantime, Eq.(201) is broken due to gluon corrections. Indeed, as a 
consequence of (|2TD,(|T2"D one has, 

Tis = -^^-(g 2 + 4m 2 )/ G (g 2 )-g 2 H(g 2 ) + o(-^), 

up = -^-g 2 (W) + s(g 2 )) + o(i). (33) 

m A z 
The axial Ward identity (for £ = 0) allows to find ft^ (Q), 
Am 2 Tlp{Q) = -Am < qq > +Q 4 Tlf(Q); 

n? } (Q)= -^(/ G (Q 2 ) + H(g 2 )). (34) 

Therefore the discrepance in eq.(201) of [8] is solely due to non-zero S, 

U P (Q) + Q 2 tlf(Q) - Tls(Q) = -4m 2 5(Q 2 ) + 0. (35) 
In the soft-momentum expansion, (£ = 0), one obtains, 

(0)(n\ N ° I 1 ( 2 A A 2 ,\ , < G lu>\ ( 1 <Gl„> 



Q 2 \ \ m? J 24N c m 2 J \Q 120iV c m 4 / 
4m 2 n P (Q) = -Am<qq> +Q 4 Tlf(Q); 

U P (Q) + Q 2 Tlf(Q)-Tl s (Q) = ^i^ + O(%)^0. (36) 

Thus these identities prove that the results of [8] should be revised. 

6. Determination of meson masses 

The invalidity of the cx-model relations is of course a bad news since we cannot 
calculate reliably two-point correlators in the GNJL model for finite momenta ~ 
lGeV but would like to obtain the information about spectra from the soft-momentum 
expansion. However, it is not a property of the Chiral Symmetry Breaking to provide 
the linear cr-model relations. The QCD-motivated pattern of the CSB [7] includes the 
radial excitations of pions and scalars with unequal masses and coupling constants 
that is an extra source of the CSB. Two interpolation schemes of mass spectrum 
determination can be developed. 
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1. Following the planar limit of the QCD, eq.([l]) one can make the two-resonance 
ansatz for scalar and pseudoscalar correlators in the GNJL model, eq. (|l9|) , 

u s (o) = < ^ >2 r- 1 i <gg> _ z q | z ° i c°- 

W) m 2 a + m Q 2 + Ml Q 2 + M 2 + ' 

T\ P (Q) - < ^ >2 F-i i <Ciq> - ^ i Z ? I C« (37) 

We remark that for the NJL-type models the constants C]_ ' =0. ^From the require- 
ment of asymptotic CS restoration (^) it follows that, 

C a = C* = C = <m> < 0; (38) 

m 

Z° + ZZ = Z% + Z1- (39) 

Z°M 2 + Z{Ml, = Z^M 2 ,. (40) 

The first two relations can be fulfilled in the conventional NJL model which corre- 
sponds to the one-resonance ansatz, Zi' n = 0, whereas the last one can be saturated 
only in a two-resonance model including at least 7r'-meson. 

The soft-momentum limit of the correlators is connected to the structural con- 
stants of the chiral lagrangian [8, 13], 

§2 + C = 8BZ(-2L B + H 2 ). (41) 

7r' 

When eliminating the unobservable constant H 2 one comes to the large- N c sum rule 
for meson parameters based on the phenomenological constant L 8 [13], 

<ya <7cr <y-K 

% + k~kr 32B '^ (42) 

Let us interpolate now the soft-momentum expansion (|20|) for r CT)7r of the GNJL model, 
eq.([5TD by means of the two-resonance ansatz (|37p, 



C 2 M 2 M 2 , Z°M±, + Z°Mt 

" Z%M 2 +Z?M*' 1 ZgM*+ZZM*' 

< = C 2 — ; al = -C 2 ^— (43) 

As a consequence, seven meson parameters M a yy, Zq'J are determined by entries 
of the GNJL models, A, m, < G 2 ^ >. The preliminary estimates show that the gluon 
condensate of the GNJL model should be taken less than the similar nonperturba- 
tive parameter in the high-energy OPE expansion [11,12]. It is not surprising since 
partially the nonperturbative gluon effects are taken into account in the four-quark 
interaction. 
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2. In the recent literature there appeared few conjectures [4,8,10] of the validity 
of (j-model relations like I n = I a (see eq. fll4]) ) and m a = 2m dyn in the case of a 
NJL model coupled with external vector fields. In particular, in the review [10] the 
influence of external electric or magnetic fields with constant strength and direction 
was considered for the case of two-point correlators. These configurations seem to 
form the same background condensate as gluons at the lowest order in fields and 
external momenta. However it can be checked that the restorating of the Lorentz 
symmetry of the vector-field condensate (by averaging over Lorentz rotations) mixes 
the longitudinal and transversal (to the external momentum) components. It results 
in the breaking of linear a-model relations (|i~4]). 

Nevertheless one can develop the soft-momentum scheme for the determination of 
scalar meson masses based on the approximation, 

T a ~ g 2 (/ G (0) + ~(0))+4m%(0); I\ ~ q 2 (J G (0) + 5(0)) ; 
, 2 4m 2 / G (0) al 

im ° = 1MTW) = ^ (44) 



that should be compared with adiabatic cr-model relations (|22|). The Nambu relation 
is not valid in this scheme, m 2 > 4m 2 . 

If to append the two-resonance ansatz to the adiabatic approximation one arrives 
to the remarkable condition for scalar meson masses, 

< = (45) 

which guaranties the validity of above approximation. As Zq > Zq one expects that 
in the chiral limit the pion and a-meson decay constants are different, F n < F a . 
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